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I. INTRODUCTION 



With the commissioning of new powerful facihties with high production rates, low- and 
medium-energy meson physics will experience renewed interest. B physics at higher energies 
and K physics at lower energies have stimulated large efforts on the experimental as well 
as the theoretical side; the field of rj and rj' physics provides a considerable amount of open 
questions, and the new facilities are expected to address them from the experimental side. 
The anticipated numbers of 10^ - 10^ observed etas per year at CELSIUS (~ 2.2- 10^), ITEP 
(~ (0.27-2.7) • 10^) and DA$NE (~ 3.2-10^) |T| will allow for experiments which on the one 
hand supply precise figures on the more frequent eta decays and which on the other hand 
focus on rare eta decays. Such rare decays can supply valuable information on anomalous 
processes or a possible C violation in eta decays. Another interesting question is whether it 
will be feasible to observe some of the rare eta decays at the new laser backscattering facility 
GRAAL at Grenoble (~ 10 decay events per second) or even at the cw electron facilities 
at Mainz (MAMI), Bonn (ELSA) and Newport News (CEBAF), where many eta photo- 
and electroproduction experiments are scheduled or already being carried through. On the 
theoretical side the development of chiral perturbation theory [0,^ as an effective theory 
for the confinement phase of QCD has supplied a consistent framework for the calculation 
of low-energy processes in the SU{2) as well as the SU{3) flavor sector of QCD. Chiral 
perturbation theory has turned out to be a valuable tool in the investigation of meson 
interactions at low energies. However, the theory seems to work much better in the SU{2) 
sector than in the SU{3) sector due to the comparatively large mass of the strange quark. 
In particular, processes involving the eta meson have confronted chiral perturbation theory 
with various problems which could only be solved partly by considering next-to-leading-order 
terms and electromagnetic corrections or by going even beyond next-to-leading order (see, 
e.g., @J^). These problems are also related to the large r] - rj' mixing angle and the t/yi(l) 
problem. In pure chiral perturbation theory up to 0{p^), the eta singlet field is integrated 
out 0, because the mass of the corresponding physical r]', rrir^i = 957.7 MeV, is larger than 
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the low-energy scale Acpt ~ fTT'p- On the other hand, current-algebra-like calculations with 
phenomenologically determined decay constants and rj - rj' mixing, which keep the eta singlet 
field as an explicit degree of freedom, yield reasonable results. Since we will restrict our 
calculation of the process 1] Tm^'-f to tree level, we will choose such a phenomenologically 
inspired approach. The charged mode of this process is of particular interest, because it can 
supply information on a Wess-Zumino-Witten contact term [^,0 involving the interaction of 
three mesons and two photons (for a review of anomalous processes see, e.g., [§]). The eta 
meson is the only particle of the pseudoscalar octet which decays through such a mechanism. 
Another possibility of testing this special kind of contact term in a scattering experiment 
is the process 77 tt+tt^tt^ which has recently been treated in chiral perturbation theory 
0. Due to the lack of meson collision facilities, it is difficult to extract information on this 
vertex from a scattering experiment with a meson in the initial state, because in practice 
such experiments involve virtual particles and necessitate uncertain extrapolations. On the 
experimental side a data analysis dating back to the sixties with rather low statistics gives 
upper bounds on the branching ratios of the charged mode 77 —* 7r+7r^77 |T^-|T^: 

[11] : r (77 ^ 7r+7r-77) /F^oK 2.1 x 10-^ (1) 
H : r (r^ ^ 7r+7r-77) /F^K 3.7 x 10"^ (2) 



In the analyses [|TT| and |]T2l the upper limits for the branching ratios were derived for a 



missing mass of neutral particles larger than 195 MeV. Hopefully, the situation will improve 
when future experiments will be carried through. 



II. KINEMATICS AND OBSERVABLES 

The four-momenta and polarization vectors for the charged decay mode rj 7r+7r~77 
are defined in Fig. 1. For the neutral decay mode we use analogous descriptors. The full 
kinematics of a decay process with four particles in the final state requires five independent 
kinematical variables (see, e.g., |TB|,|T^). For the definition of these variables we will consider 
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three reference frames: the rest system of the eta meson E^, the dipion center-of-mass system 
E^TT) s-nd the diphoton center-of-mass system E^^. Our kinematical variables are (see Fig. 2) 

• Stt, the square of the center-of-mass energy of the pions, 

• s-y, the square of the center-of-mass energy of the photons, 

• the angle of the pion with momentum ki in Etttt with respect to the direction of 
flight of the dipion in E^, 

• , the angle of the photon with momentum qi in E-y^ with respect to the direction 
of flight of the diphoton in S^, 

• 4>, the angle between the plane formed by the pions in E^ and the corresponding plane 
formed by the photons. 

In order to define these variables more precisely we introduce a unit vector v along the 
direction of flight of the dipion in E^, and unit vectors c and d along the projections of ki 
perpendicular to v and of perpendicular to —v, respectively, 

1/2 



= (^ki — vv ■ kij I ki — {k,\ ■ i?j 



(3) 



-1 1/2 

d^{qi-vv- qi) / [q^ - (g! • vf\ . (4) 
With these definitions the five kinematical variables are defined as follows: 

Stt = {h + hf , (5) 

S7 = (?i + ^2)^ , (6) 

cos 6*^^ = V ■ ki/ \ki\, (7) 

cos 61^1 ^ -V ■ qi/ \qi\, (8) 

cos 4> = c ■ d. (9) 

The physical region of the decay process is refiected in the range of the kinematical variables: 
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< s^<{mr,- 2m^f , (10) 

4m^ < s^< (m^ - ^) , (11) 

0<^.i<7r, (12) 

0<^-yi<7r, (13) 

< < 27r. (14) 

The invariant matrix element squared, | will be expressed in terms of Lorentz scalar 

products of the five momenta ki,k2,qi, q2 and p. One of these momentum vectors can 
be eliminated because of momentum conservation. In order to express the Lorentz scalar 
products in terms of the kinematical variables specified above we now introduce adequate 
linear combinations of the momenta: 

K^ki + k2, (15) 

L = ki- k2, (16) 

Q^qi + q2, (17) 

R = qi- q2. (18) 

For further reference we need the expressions 



K- 


K = St,, 


(19) 


Q 




(20) 


K 




(21) 


K 


R = X cos 9^^ , 


(22) 


L 


Q — a-j^x cos 6t,^ , 


(23) 


L 


■ R = aT,[K ■ Q cos 0,,^ cos 





{st^s^Y^'^ sin sin 9^^ cos , (24) 



with 



X = ^J{K-Qf -s^s^, (25) 
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a, = ^Jl - Aml/s,. (26) 
The ten Lorentz scalar products in \M.\'^ can now be expressed as follows: 





k2 = 


\(K-K-2ml), 




(27) 




■ q2 = 


\q-q. 




(28) 


h 


• Qi = 


^{K-Q + L-Q + K-R + L- 


R), 


(29) 


k2 


• Qi = 


^{K-Q-L-Q + K-R-L- 


R), 


(30) 


ki ' 


■ Q2 = 


-iK-Q + L-Q-K-R-L- 


■R), 


(31) 


k2 


■ ?2 = 


^{K-Q-L-Q-K-R + L' 


R), 


(32) 


P- 


k^^ 


ki ■ qi + ki ■ q2 + ki ■ k2 + ml, 




(33) 


P- 


k2 = 


k2 ■ qi + k2 ■ q2 + ki ■ k2 + ml, 




(34) 


P 


■ 91 = 


qi-q2 + ki ■qi + k2- qi. 




(35) 


P 


• ?2 = 


qi- q2 + ki- q2 + k2- ?2- 




(36) 



The differential decay rate can be written as 

d^r {rj 7r7r77) = 2-^\-^m-^C-^a^x \M\^ 

ds^ds^dcos decs 9^-^ dcf), (37) 

where the symmetry factor C is equal to 4 in the decay rj —>■ tt'^tv^jj because of two pairs 
of identical particles in the final state, and equal to 2 in the decay rj 7r"'"7r~77 because of 
one identical particle pair in the final state. Now we will proceed to investigate how chiral 
dynamics manifests itself in the Lorentz-invariant matrix element Ai. 

III. CHIRAL DYNAMICS OF THE DECAY r? ^ tttt-/-/ 

We will restrict our calculation of A4 to the leading-order contributions. Since the 
process rj — > 777777 involves the electromagnetic interaction of an odd number of pseudoscalar 
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mesons, the leading contributions must contain a vertex of odd intrinsic parity. Such a 
vertex is at least of O(p^) in the momentum expansion, and thus, according to Weinberg's 
power counting [Q], we expect the leading contribution to be of O^p"^). Consequently, the 
interaction Lagrangian we will use for our tree-level calculation contains the standard 0{p^) 
piece and the anomalous Wess-Zumino-Witten Lagrangian 1^,0], but no terms from the 
Gasser-Leutwyler Lagrangian of 0{p'^): 

L- — L- -h f-wZW 

= ^tr{{D'^UyD^,U) + ^tr(x^f/ + xU^) 
eN 

ie'^N 1 1 

''-e^"'''^d^AMT{Q^iUdfsU^ + dpU^U) - -QU^QdpU + -QUQdpU^), (38) 



.(4) 



where £0123 = 1- In Eq. ( pSD we have only listed those terms of jCwzw which actually give 
a contribution to the invariant amplitude. The covariant derivative is defined as 

D^U = d^U + leA^ [g, U\ , (39) 

where the matrix Q represents the electromagnetic charges of the three flavors in S'f/(3), 

Q = diag(2,-l,-l)/3. (40) 

The matrix 



X = 2Bom (41) 



contains the quark masses. 



m = diag(m„,md,mj , (42) 



where Bq is related to the quark condensate and is given by the relation (m^ + rrid) Bq = m^. 
The meson field operators are represented by the matrix U = exp (i$/F$), where the nonet 
field matrix $ can be decomposed into an octet and a singlet part, $ = + with 
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and 



V ->y 



$1 = W-r/odiag (1,1,1) 



(43) 



(44) 



In the expansion of U the decay constants F^, Fg or Fq will be inserted for the constant 
depending on whether the constant belongs to a vr, r/g or rjo field. We will use = 
93MeV, Fg = 1.25F,r and Fq = I.OGFtt as numerical values [|15[. Our calculation in chiral 



perturbation theory will be carried out with the group theoretical octet and singlet eta 
states, |?7g) and |?7o). We will introduce rj-i]' mixing via the phenomenological mixing angle 



e = -20° 0: 



|r^) = cos6'|r7g) - sm9\r]o), 
\ri') = sin6'|?78) + cos6'|r7o). 



(45) 
(46) 



The Feynman diagrams contributing to the charged decay mode i] 7r+7r^77 of 0{p^) are 
displayed in Fig. ^ There are three different classes of diagrams at tree level: diagrams with 
a four-meson vertex of O(p^), a propagating neutral meson, and a decay vertex into two 
photons of 0(|3^) (class 1), Wess-Zumino-Witten contact terms of 0{p'^) (class 2) and internal 
bremsstrahlung diagrams, where one photon is emitted off a charged pion line (classes 3.1 
and 3.2). The first class of diagrams is gauge invariant by itself, whereas the amplitudes 
corresponding to the diagrams from the second and third class have to be added in order 
to obtain gauge invariance. In the neutral decay mode r] 7r°7r°77 only the first class of 
diagrams is relevant (Fig. 

Starting from the general chiral Lagrangian of Eq. (|38|) , we now list the interaction terms 
relevant for the process rj 7r+7r^77: 



(4),2730 
WZW 



12v^7r2F2 



X (^dpr^s + ^ dpvo] , (47) 



Fs ^ '° F 
J-wzw - i2v/37r2F2 ' 



X (^dpr]s + ^ 9;3r/ol , (48) 



^ '° F 



^(2),4<A ^ -^^^-Sq ( 4 (m^ + mrf^ 
24 1, F2F| 



8a/2 (m„ + m^) + _ 
8 (m„ + md) + _ , 
8 (m„ - m^) + _ 

-r^gTT TT TT 



8^2 (m„-md) + _ o\ /cn^ 
H ^ — = rynvr^TT TT , (50 

£(2),i72</. = ieA" (dp'K-'K+ - dy-R-) . (51) 

The invariant amplitude for rj 7r+7r~77 is then given as the sum of the amphtudes from 
the three classes of Feynman diagrams (Fig. 

M=Mi + M2 + M3, (52) 

where 

ie'^Nc ^^^0 I Bq {m^ + rrid) /cos 6* y2sin6'\ 



^ 45o (m„ + mrf) f cos 9 ^/2sm9' 



3{2q,.q2-ml)\FlFs F^ ) 
Boimu — rriii) /cos 6* ^/2sm9^ 



3 {2q, ■ q, - ml^) \F^Fs F^Fo 



(53) 



_ je^iV^ . «^ /^"^^'^ y2sin^^ 
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ie^Nc auaB ( ^^^^ sin 6' \ [ £2-^1, , Si ■ ki 

+£i,fiqi,aki,uP(3 r- + £2,fMq2,aki,uP/3 r- > ■ (55) 

?2 ■1^2 qi ■ 1^2 ) 

The masses in the propagators can be expressed in terms of the physical masses using the 
relations 



m. 



m. 



vo 



cos^ 9 + m^, sin^ 9, (56) 
m.^ sin^ 9 + m^, cos^ 9, (57) 



ml, = mlo. (58) 

As we will not be concerned with photon polarizations in the final state, we carry out the 
sum over the polarizations of the real photons, 

\M\''= M{K^,K2)M*iKi,K2). (59) 

For that purpose we exploit the completeness relation 

^7 ^5 (/«) -9js, (60) 

K 

which is based on current conservation (see, e.g., [|T6| ): 



M,{k,,K2)M*{k,,K2) = Y -^^-^fXA(«:iXp(«:i)^2,.(«:2) 4.(^:2) 

Kl,K2 

gx,g„rM':''MY* = M^^Mx./ e {1, 2, 3}) . (61) 



,Xa \Apr* _ A^Acr 



The result for | p is too complicated to be displayed in this contribution. It depends on 
all Lorentz scalar products which can be constructed from the four final-state momentum 
vectors ki, k2, qi and q2- 

For the neutral decay mode 77 — > 7r'^7r'^77 we obtain the interactions 

. -I^,e^-%A^A^ [^^Oy + + . (62) 
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+ F-F.r, """" " + V3F|Fo " " 

+i^%i^*%-°.° + i^^^^»*,„.„). (63) 

which lead to the diagrams in Fig. ^ and resuh in the Lorentz-invariant matrix element 



^ 8 (m„ + m^i) ( cosO y^sin^' 



32 (m„ + m^) / 008 6* \/2sm.O' 



3v^(2gi-g2-m2 ) VF^Fq fI 



(64) 



Summing over the possible photon polarizations ni and ^2 we obtain a compact result for 
the invariant matrix element squared, 

,-^,9 e^N^B?. , .0 I rriu — rrid ( cos9 v^sin^^' 



TUu + 171(1 f cos9 V2sm9\ 



3^/3 (2g 



2 U F| FiFn 



1 ■ g2 - \ -^8 / 

4 (rriu + m(i) / cos 9 a/2 sin 6* ' 



2 



3v^(2gi-g2-m2j VFsFo F^^ , 
We note that the final result only depends on one Lorentz scalar product, namely qi ■ q2 



(65) 



IV. RESULTS AND DISCUSSION 



Having determined we proceed to investigate the decay spectra by integrating Eq. 

(P7|) numerically. First we will discuss the Dalitz plot d'^T/ds-^ds^ for the charged decay 
mode. We note that the soft-photon limit for any of the two photons implies 0. 
The contour plot for the full chiral perturbation theory calculation (Fig. |^) clearly shows 
the infrared bremsstrahlung singularity for ^ 0. The Dalitz plot also refiects a pole 
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at s-y = m^o = 18219 MeV^ due to the class 1 Feynman diagram with a propagating 7r° 
(see Fig. ^). In the neighborhood of this pole it is impossible to distinguish between the 
processes rj — > 7r+7r^77 and t] tt+tt^tt^. If we switch off the diagrams of class 1, the 
invariant amplitude is still gauge invariant, and the corresponding Dalitz plot (Fig. ^ is 
very similar to the full calculation (Fig. |^) except for the region around the 7r° pole. We 
conclude that the diagrams of class 1 do not contribute significantly to the process which 
also becomes evident from the diphoton spectrum dT/dz (Fig. |^), where z = s^/m^. The 
effect of the 7r° pole is confined to a very small region around = m'^o- The calculation 
without class 1 diagrams almost coincides with the full calculation except for this region. 
The reaction mechanisms of class 2 and class 3 diagrams dominate the spectrum over a 
wide energy range. We conclude that the detection of this decay mode should be a good 
indication for the presence of a Wess-Zumino-Witten contact term (class 2). Towards small 
values of s-y the bremsstrahlung diagrams (class 3) will be responsible for a divergence in the 
spectrum. A measurement of the steep rise in the spectrum for vanishing depends on the 
resolution of the detector facilities or, more accurately speaking, on the minimum photon 
energy detectable. Hence, in both experiment and theory it is only possible to determine 
the partial decay rate as a function of an energy cut 6mbrems applied around y/s^ = 0. We 
obtained T [t] —>■ 7r^7r^77) by integrating the diphoton spectrum for the calculation without 
class 1 diagrams. According to Fig. ^ the error introduced by this approximation should be 
negligible. In Fig. || we show the result as a function of Smtrems- Comparing our absolute 



numbers for the partial decay rate with the total eta decay rate Ttot = 1-2 x 10 ^ MeV [ITO 



we note that the resulting branching ratio is well within the reach of the facilities mentioned 
in the introduction. Using an energy cut at Srribrems = 195 MeV, our result for the partial 
decay rate is smaller than the experimental upper limits in Eqs. (|l]) and (0) by about three 
orders of magnitude. 

Let us now turn to the neutral decay mode rj iT^TT^'-f'-f. In contrast to the charged decay 
mode, where the diagrams of class 1 can be neglected, these diagrams generate the only tree- 
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level contributions in the case of neutral mesons. For this reason we want to analyse the 
structure of the vertices in the diagrams more closely. Whereas the two-photon-one-meson 
vertices have been investigated extensively in the two-photon decays of the ir^, rj and rj', the 
four-meson vertices, in particular the rjsrjsn^Tr'^ interaction, are not known with comparable 
precision. However, these vertices are of great interest, because they are directly related to 
the sum or difference of the light quark masses (see Eq. (|63|) ). While the r^vrvrvr and ri'Tjirir 
interaction can be directly investigated in the decays rj tttttt and t]' rjTnT, this is not 
the case for the Tjijirir interaction. Moreover, the minimum center-of-mass energy for the 
scattering process vrvr —>■ rjr] is beyond the convergence radius of chiral perturbation theory, 
and TTr] tti] scattering is not a realistic alternative either. As a consequence, the only 
possibility to get information on the rjrjiT'K interaction is the investigation of a composite 
process with at least one additional vertex. At first glance the decay rj tt^tc^'^'j seems 
to be a good candidate since G parity is conserved at the r^r^vr^vr^ vertex, whereas in the 
diagram with a propagating the corresponding vertex violates G parity and vanishes in 
the isospin limit. However, as is well-known from the decay r] tttttt, G parity is violated 
significantly and isospin breaking is reflected by the fact that the rjinnT vertex function is 
directly proportional to the quark mass difference (see Eqs. ( [5D| ) and (§^)). For this reason 
the diagram with the propagating neutral pion cannot be neglected. Since m^o lies within 
the range of integration for s^, the diagram with a propagating tt*^ will cause a pole in 
the amplitude M.. Thus, despite G parity and isospin violation, this diagram is strongly 
enhanced in comparison with the diagrams with a propagating 773 and tiq. The diphoton 
spectrum clearly demonstrates the fact that the rjsrjsiT^TT^ interaction plays a minor role in 
the process (see Fig. |^). 

In the pole region it is impossible to distinguish the decay mode r] tt°tt'^77 from 
7] TT°TT°TT°. Thus, in order to obtain new information on chiral dynamics from the decay 
1] tt°tt'^77, one should choose an energy regime for which is sufficiently far away from 
this pole. 

At 0{p'^) effective chiral Lagrangians do not reproduce the experimentally determined 
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interaction strength of the r^vr'^vr^vr^ and rj'rjTi^TT'^ vertices (for the correct treatment be- 
yond tree level see Therefore, we have also performed a phenomenological calculation, 
where we fixed these interactions using the experimentally observed branching ratios 
T{ri^ 7r07r°7r°) /Ttot = 0.319 and T [r]' rjn^TT^) /Ttot = 0.208. Whereas experimental data 
agree very well with a constant tree-level prediction for rj ti^ti^tt^, such an assumption 
does not seem to work so well for t]' rjir^Tc^ (see, e.g., the discussion in ^^). However, 
this is of only minor importance in our case, because the pole at m^, in the amplitude is 
far outside the physical region of Sj. It turns out that the result with the vertices fixed by 
experiment is larger than the prediction of chiral perturbation theory by about one order of 



magnitude. Furthermore, we also found that n-i] and n-i]' mixing according to [|l^ is neg- 



ligible in this process. Future work will focus on higher-order corrections to our tree-level 
calculation. 

The question whether it is realistic to measure the decay r] — > 7r°7r°77 must be decided 
from the partial decay rate T [r] ^ 7r°7r'^77), which we obtain by integration of the diphoton 
spectrum applying a cut of the width 26m around the pole at m^o (Fig. From this 
figure it becomes obvious that this partial decay rate is rather small, but possibly within 
the reach of the new eta facilities. However, it will probably be impossible to draw any 
precise information on the G-parity-conserving Vj^rfTX^Ti^ coupling from measuring this pro- 
cess, because the corresponding Feynman diagram gives only a small contribution to the 



complete amplitude. Finally, we want to mention that the calculation in ||T9[ based on old 
experimental data yields F (77 — > 7r'^7r'^77) = 1.38 x 10^^ eV. 

We conclude that the rare eta decay r] 7r+7r^77 is an interesting test of the anomalous 
Lagrangian, because it is a simple way to access the three-meson-two-photon vertex predicted 
by and ||^ and allows for a consistency check with the results from 77 tt+tt^tt^. The 
neutral decay mode also investigated in this contribution will be much harder to detect. 
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FIGURES 



7i {Qi,£i (ki)) 




FIG. 1. Notation for the kinematics of the process ij 7r"'"7r 77. 
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FIG. 2. Choice of kinematical variables in the four-particle decay rj tttt^^. 
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Class 3.1 Class 3.2 

FIG. 3. Feynman diagrams of the process tj tt'^tt^j^. The numbers in the interaction blobs 
denote the order of the vertex in the momentum and quark mass expansion. 
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FIG. 4. Feynman diagrams of the process rj tt tt 77. 
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FIG. 5. Dalitz contour plot for the full calculation in chiral perturbation theory. The 
distance between two contour lines corresponds to a difference of 1.47 x 10~^^MeV~^ in the 
doubly differential decay rate d'^T/dsT^-dsj. The outermost contour line denotes the level 
d^r/ds^ds^ = 1.467 X 10-^'^MeV-^. The dotted area is the physical region of the process. The 
dots in the s-j^-s^ plane represent the points where the invariant amplitude was evaluated. 
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FIG. 6. Dalitz contour plot for the chiral perturbation theory calculation without class 1 pole 
diagrams. The contour spacing is the same as in Fig. ||. 
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FIG. 7. Diphoton energy spectrum dT/dz {z = s-y/m^) for the decay 77 — 7r"'"7r'~77. The 
dotted line is the full calculation, the dashed line is calculated without class 1 diagrams and the 
solid line is proportional to the phase space integral. 
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FIG. 8. Partial decay rate T{ri — > tt'^tt 77) as a function of the energy cut Smbrems around 
the bremsstrahlung singularity at = 0. 
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FIG. 9. Diphoton spectrum for the decay rj — 7r°7r°77 {z = s^/rn^). The long dash-dotted 
line is the prediction of chiral perturbation theory, the short dash-dotted line is the calculation 
in chiral perturbation theory, where the n^n^ijs'ijs interaction is set to 0, and the dotted line is 
the prediction of chiral perturbation theory for m„ = m^. The short dashed line is a calculation, 
where we determined the strength of the r77r^7r*'7r^ interaction from the decay ry tt^tt^tt^ and the 
■qri'iT^'K^ interaction from the decay Vj' — > t/tt'^tt^. In the long dashed curve an additional rygr/sTr'^Tr^ 
interaction is included. The solid line is proportional to the phase space integral. 
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FIG. 10. Partial decay rate for the decay rj 7r°7r°77 as a function of the energy cut Sm in 
Sj^^ around m^o. The long dashed line is the calculation with the r/Tr^Tr'^Tr^ and vi'-qir^-K^ vertices 
fixed by experiment. For the solid line we added a rj^rj^TT^TT^ interaction from chiral perturbation 
theory. 
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